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INTRODUCTION 

Late-type stars radiate most of their energy 
in the near- and mid-infrared regions. The en- 
ergy distribution of Mira variables peaks 
around 2 /t m and the well-known infrared 
source, IRC + 10216, is very bright between 
2 and 20 /tin. Many infrared stars similar to 
IRC + 10216 are believed to be long-period 
variables, their thick circumstellar shell almost 
totally obscuring the central star. Reemission 
of stellar radiation by warm dust grains causes 
the infrared continuum flux. The extended 
shells of both visible and unidentified infrared 
cool stars also emit molecular emission lines at 
centimeter and millimeter wavelengths. Since 
radio lines can be excited by infrared photons, 
a combination of infrared and radio observa- 
tions is very useful in determining molecular ex- 
citation mechanisms (Kwan and Scoville, 1974; 
Deguchi and Iguchi, 1976; Elitzur et al., 1976; 
Bujarrabal et al., 1980). 

Late-type stars are characterized by mass 
loss (see Goldberg, this volume). Sporadic ejec- 
tion of matter or modulation of a continuous 
process may produce stratification of the cir- 
cumstellar shell (Bernat, 1981; Ridgway, 1981). 
As a result, molecular lines can serve to probe 
the physical conditions in different layers. In 
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particular, SiO masers (rotation lines in the 
ground and excited vibration states) and in- 
frared vibration/rotation molecular lines, 
which need extreme excitation conditions (i.e., 
high gas densities and temperatures), occur 
close to the stellar photosphere. On the other 
hand, millimeter thermal emission from CO 
and linear carbon chain molecules — the cyano- 
polyyne family (HC 2n + ,N)— takes place in the 
circumstellar envelope at several 10 to 10 3 
stellar radii. Different shell layers can therefore 
be sampled by observing appropriate molecular 
transitions (Figure 3-1). 



Figure 3-1. Schematic drawing of the 
circumstellar envelope. 
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CIRCUMSTELLAR MOLECULES 

Molecules in space were discovered through 
the ultraviolet absorption lines of CH, CN, and 
CH + in stellar spectra (Swings and Rosenfeld, 
1937; McKellar, 1940; Adams, 1941). The first 
detection of interstellar molecules at radio 
wavelength, the OH line at 18 cm, was made 
by Weinreb et al. in 1963. The improvement of 
radio astronomical techniques resulted later in 
the detection of polyatomic molecules, NH 3 
and H 2 0, around 1.3 cm (Cheung et al., 1968, 
1969). So far, about sixty molecules and many 
of their isotopes have been detected, mostly in 
the millimeter range, in the star-forming region 
of Orion, in the galactic center (Sgr B2), and 
in the circumstellar envelope of the carbon star 
IRC + 10216, as well as in a number of dark 
clouds. 

The first discoveries of circumstellar 
molecules were made through the detection of 
OH and H 2 0 maser lines in oxygen-rich late- 
type stars, mainly Mira variables and red 
supergiants (Wilson and Barrett, 1972; Knowles 
et al., 1969; Schwartz and Barrett, 1970; Robin- 
son et al., 1971; Nguyen-Q-Rieu et al., 1971). 
SiO maser emission was later detected in these 
objects (Snyder and Buhl, 1974; Kaifu et al., 
1975; Spencer et al., 1981). Subsequent sys- 
tematic surveys in the 1612-MHz OH line in the 
galactic plane resulted in the detection of OH 
sources whose spectra exhibited the charac- 
teristics of OH stellar masers (Winnberg et al., 
1973; Johansson et al., 1977; Bowers, 1978; 
Baud et al., 1979a, 1979b). These OH sources, 
which also have infrared but no visible counter- 
parts (Schultz et al., 1976; Evans and Beckwith, 
1977; Glass, 1978; Epchtein and Nguyen-Q- 
Rieu, 1982), form the so-called “unidentified 
OH/IR” class of objects. About 300 stars and 
unidentified IR objects are known to emit at 
least one of these maser lines (Engels, 1979). 

Thermal CO emission has been detected in 
the envelope of carbon-rich as well as oxygen- 
rich stars (Solomon et al., 1971; Zuckerman et 
al., 1977). Other molecules containing carbon, 
such as HCN, HNC, and HC 3 N have also 
been detected in some carbon-rich stars (Jewell 


and Synder, 1984; Olofsson et al., 1982b). In 
particular, 22 molecular species and their iso- 
topes have been found in the envelope of IRC 
+ 10216, including a complex member of the 
cyanopolyyne family, HC U N, which is also 
the heaviest molecule so far detected in our 
Galaxy (Bell et al., 1982). Recently, Thaddeus 
et al. (1984) have identified nine of the pre- 
viously unassigned lines in IRC + 10216 as 
those emitted by the SiC 2 radical, which turns 
out to be a compact symmetric ring molecule. 
Table 3-1 lists the circumstellar molecules. (See 
also Glassgold and Huggins, this volume.) 

Circumstellar molecular lines are char- 
acterized by their broad line width, which 
reflects the large-scale expansion of the shell 
(-5 to 50 km s' 1 ). 

MASER EMISSION 

The radio emission from OH, H 2 0, and 
SiO presents characteristics of nonthermal pro- 
cesses. The brightness temperature of maser 
sources can be as high as 10 14 K. While OH 
and H 2 0 masers correspond to transitions in 
the ground vibration state v = 0, SiO masers 
occur in excited vibration states u = 1 , 2, and 
3. Recently, SiS maser emission in the ground 
state has also been detected in the envelope of 
IRC + 10216 (Henkel et al., 1983). 

Masers are variable sources whose intensities 
are correlated with those observed in the near 
infrared (Hjalmarson and Olofsson, 1979). A 
high degree of polarization has been observed. 
Population inversion can be achieved through 
pumping by infrared radiation from warm cir- 
cumstellar dust grains (Elitzur et al., 1976; Bu- 
jarrabal et al., 1980) or from the central star 
(Kwan and Scoville, 1974; Deguchi and Iguchi, 
1976; Bujarrabal and Nguyen-Q-Rieu, 1981). 
Details of the pumping mechanisms will be 
discussed in the following sections. 

Two-Level Maser 

It is possible to describe a simplified maser 
theory by representing all details of the pump- 
ing mechanism in a two-level maser using the 
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Table 3-1 
Molecules in Stars 


Molecule 

Transition 

Star 

Characteristics 

OH 

Radio 

Oxygen-rich 

Maser emission 

H 2° 

Radio 

Oxygen-rich 

Maser emission 

SiO 

Radio 

Oxygen- and carbon-rich 

Maser and thermal 
emissions 

CO 

Radio 

Oxygen- and carbon-rich 

Thermal emission 

SiS 

Radio 

Carbon-rich 

Maser and thermal 
emissions 

CN 

Radio 

Carbon-rich 

Thermal emission 

cs 

Radio 

Carbon-rich 

Thermal emission 

HCN 

Radio 

Carbon-rich 

Thermal emission 

HNC 

Radio 

Carbon-rich 

Thermal emission 

HC N 

Radio 

Carbon-rich 

Thermal emission 

hc 5 n 

Radio 

Carbon-rich 

Thermal emission 

hc 7 n 

Radio 

Carbon-rich 

Thermal emission 

hc„n 

Radio 

Carbon-rich 

Thermal emission 

C 2 H 

Radio 

Carbon-rich 

Thermal emission 

C 4 H 

Radio 

Carbon-rich 

Thermal emission 

C 3 N 

Radio 

Carbon-rich 

Thermal emission 

ch.cn 

Radio 

Carbon-rich 

Thermal emission 

SiC 2 

Radio 

Carbon-rich 

Thermal emission 

nh 3 

Radio and IR 

Carbon-rich 

Thermal emission 

C 2 H 2 

IR 

Carbon-rich 

Thermal emission 

C 2 H 4 

IR 

Carbon-rich 

Thermal emission 

CH 4 

IR 

Carbon-rich 

Thermal emission 


pump rates (Litvak, 1972). The populations of 
the upper and lower levels n and n ( , are given 
by the statistical equilibrium equation: 

+ B UI + C », + P J 

= n l (B/ u OI/(4r) + C lu + PJ , 

where A „ B , are the Einstein coefficients for 
spontaneous and stimulated emissions, B fu is 
the absorption coefficient, and C ( and C /H are 
the collision rates. The pump rates, P and 
P /m> correspond to indirect population transfer 
between levels u and / through other higher 
levels which are not considered. P #/ and P /u are 
the pump rates which transfer population to the 
lower and to the upper level, respectively. / is 
the specific intensity, and U is the solid angle 
of the microwave emission. 


To simplify the explanation, we shall ignore 
spontaneous emission and collisions and as- 
sume that the two maser levels have the same 
statistical weight. The rate of fractional popula- 
tion inversion can then be derived from Equa- 
tion (3-1): 

n u ~ n t A P 

TTT^ “ P + 5 ; f2//(2rr) ’ (3 ' 2) 

where P = P f + P (u is the total net pump rate 
and A P = P f - P u( is the difference in pump 
rates to the upper and lower levels. 

Population inversion occurs when A P > 0. 
Furthermore, the maser is saturated when the 
stimulated emission comes very much into 
evidence, B u{ £lll2n » P. The population in- 
version then decreases inversely proportional to 
the stimulated emission rate. In the saturated 
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regime, the pumping is the most efficient 
because the population of the upper level is 
transferred as fast as possible to the lower level 
through stimulated emission. The saturation in- 
tensity is defined as: 


I s = 2aP/(B ug a) . (3-3) 


The population inversion rate A n = n u ~ n ( 
can be expressed in terms of the unsaturated 
rate A tx q , which is equal to (A P/P) ( n u + n ( ): 


An o 

1 +(///) 


(3-4) 


The maser intensity / can be obtained by the 
one-dimensional radiative transfer equation: 


= - kI + e . (3-5) 

dx 


The absorption and emission coefficients, x 
and €, are defined as: 

x = (n { - n u ) B ul hu/( 4nAv) , 

e = n u A hv/{ 4 ttAp) , 

where A is the spontaneous emission coeffi- 
cient, and hv is the level energy separation of 
the maser transition. 

The solution of Equation (3-5) is: 

I =J e -*x + _ L (i - e “**) (3-6) 

° K 


When population inversion occurs, ( n u > 
/? y ), x is negative, and the maser intensity in- 
creases exponentially with path length. In the 
saturated regime, Equation (3-5) gives: 


dl K o J 

""W 

dx~ 1 + (///,) 


(3-7) 


and 


I - I S f 1 ” (X ’ X s^ ’ 

where x o = x (1 + I/I s ) and x s is the length 
where saturation begins to occur. 

Equation (3-7) shows that, when the maser 
becomes saturated, the intensity increases 
linearly with path length. Litvak (1971) has 
shown that an unsaturated maser core can be 
formed at the center of a spherical saturated 
maser cloud. 

Multilevel Radiative Transfer 

A detailed study of molecular excitation re- 
quires the resolution of the radiative transfer 
and statistical equilibrium equations, involving 
many molecular levels. For an expanding 
spherical circumstellar cloud, the population of 
a level / at a distance r from the central star is 
given by: 

n \r) L (A . + B J + C ) = 

/ v ' j u u ij u 

(3-8) 

L. n. (r) (A.. + B J + C ) , 

j j y ’ v ji u u jt 7 


with 

A = 0 

for j 


u 

and 

A — 0 
ji 

for j 


A and B are the Einstein coefficients, and C 
represents the collision rates. The mean radia- 
tion field 7, integrated over the line profile and 
averaged over the angle, is usually calculated 
by use of the escape probability formalism (see, 
for example, Castor, 1970): 

J = 1 1 - m | S(r) + 0 c (r) /. . (3-9) 
The source function S(r) is: 



where g. and g. are the statistical weights, (3 is 
the probability that a photon emitted at r 
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escapes without being absorbed by the cloud, 
and @ c is the probability of a photon which 
escapes at r and strikes the central star. In 
general, /3 r , / is the contribution to the inten- 
sity of all continuuum sources: 



1 - e- 7( ^> 


rip) 


d M , 


(3-11) 


and 


|3 — — 
c 4 n 


AO ( 1 - e -T(M=1 * | 


t(m= 1) 


(3-12) 


where fx is the cosine of the angle 6 between a 
radius and the line of sight (Figure 3-2). Aft /47r 
is the dilution factor of the stellar radiation 
field (Bujarrabal and Nguyen-Q-Rieu, 1981). 

The opacity is: 



Circums-tellar 

envelop® 



1 +(7 - 1)M 2 


(3-13) 


where y is the logarithmic velocity gradient 
d\nV/d\nr , t q is related to the population in- 
version An by t q o c rAn/V(r) (Castor, 1970), 
and V(r) is the expansion velocity. 

The population distribution of the entire set 
of molecular levels can be calculated by sol- 
ving Equation (3-8) and using Equations (3-9) 
through (3-13). In particular, the population in- 
version of the maser transition i — j can be 
determined, and hence the optical depth, as well 
as the excitation temperature, T : 

n. g: 

— = — <r hv l K T ex . (3-14) 

n i 8 i 

The brightness temperature of the maser 
source can then be calculated by an integration 
along the line of sight: 

r, = {i-n,} (■-«•') • ' 3 I5 > 

with 


Figure 3-2 . Spherical model of an OH cir- 
cumstellar maser . The maser emission is con- 
fined within a double cone (hatched area) 
whose axis is along the line of sight (top). A 
schematic maser spectrum is drawn (bottom). 
The regions of the spectrum emanating from 
the hatched areas in the circumstellar envelope 
are also represented (hatched areas). 


T = 


hv 

k 


^hv/kT 


l) 1 . 


T f)h is the 3 K cosmic background temperature. 
The convolution of T a with the antenna beam 

D 

will give the observed main-beam brightness 
temperature. See Equation (3-16). 


Characteristics of Circumstellar Masers 


Circumstellar maser emission, namely that 
of OH, is often characterized by a double- 
peaked spectrum (Figure 3-3). This shape can 
be explained in terms of a spherical model in 
which the amplification path is maximum in the 
direction of the line of sight which intersects 
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the central star (see Nguyen-Q-Rieu et al., 
1979). In an expanding spherical shell for which 
y is constant and positive, the velocity field 
can be expressed as: K(r) = V (r/R y ), where 
V is the expansion velocity at the outer 
radius R. Equation (3-13) shows that the op- 
tical depth along the z axis ( 6 = 0) is t /y (see 
Figure 3-1). In the case of an expansion due to 
radiation pressure, 7 is small (7 <SC 1) and 
r t . Hence, the maximum amplifi- 

cation, corresponding to the two spikes in the 
maser spectrum, occurs in a narrow double 
cone whose apex is the central star and whose 
axis is the line of sight (Figure 3-2). The emis- 
sion at the line center corresponds to the part 
of the envelope outside the cone and perpendic- 
ular to the line of sight with the smallest path 
length. This geometrical effect is illustrated in 
Figure 3-2, in the case of a circumstellar shell 
expanding at a nearly constant velocity. The 
angle $ e of the cone corresponds to the direc- 
tion in which the line intensity drops to 1/e of 
the peak intensities, I Mq \ 

cos 0 ~ (V - AV )/V 

e v exp e ' exp * 

where AK is the line width at 1/e I Max of the 
narrow maser spikes. 

Figure 3-3 shows the 1612-MHz OH spec- 
trum and the maps of the maser emission from 
the circumstellar shell of the unidentified 
OH/IR object, OH26.5 + 0.6, obtained with 
the very large array (VLA) by Baud (1981). The 
maps indicate that the two maser peaks arise 
from two compact regions at the center of the 
source, while the emission toward the center of 
the spectrum is more extended and exhibits a 
ring-like structure. This observational result is 
consistent with the model of an expanding 
spherical shell described above. However, this 
simple model does not account for some details 
which can be explained by large-scale nonradial 
motion or turbulence which may exist in the 
shell. 

In the framework of a symmetrical model, 
the OH blue-shifted and red-shifted wings cor- 
respond to the front and back of the shell, 
respectively. Thus, the blue-shifted wing should 


lead the red-shifted wing to the observer. This 
results in a phase lag between the blue-shifted 
and red-shifted components. The phase lag, 
which corresponds to the light travel time across 
the entire shell, is - 1 month for a shell size -8 
x 10 16 cm. 

The determination of the phase lag, and 
hence the linear source size, by means of a 
monitoring of the OH maser emission, com- 
bined with interferometric measurements of the 
angular size, constitutes a powerful method to 
estimate the distance of the stars (Shultz, et al., 
1978; Jewell et al., 1979; Booth et al., 1981; 
Baud, 1981; Herman, 1983). This method has 
been extensively used to measure the distance 
of unidentified OH/IR sources for which no 
period/luminosity relation in the visible is 
available. 

Another interesting feature of the OH and 
SiO circumstellar masers is that the maser emis- 
sion is correlated with the infrared and optical 
phases (Harvey et al., 1974; Fillit et al., 1977; 
Jewell et al., 1979; Hjalmarson and Olofsson, 
1979). However, the SiO fluxes appear to have 
a phase lag with respect to the visual light curve 
(see Hjalmarson and Olofsson, 1979). The cor- 
relation with the infrared fluxes suggests that 
the maser emission is pumped by a radiative 
process. 

The pumping by 35-jtm photons emitted by 
circumstellar grains may be operative in the 
shell and gives rise to OH masers (Elitzur 
et al., 1976; Bujarrabal et al., 1980). It has been 
suggested that direct stellar radiation can pump 
OH masers (Cimerman and Scoville, 1980). The 
pumping of the SiO and SiS masers is achieved 
by the infrared photons emitted by the central 
star. Although the OH maser emission extends 
far from the star, at several 10 3 stellar radii 
(Reid et al., 1981), SiO masers occur very close 
to the stellar atmosphere at a few stellar radii 
(Elitzur, 1980; Bujarrabal and Nguyen-Q-Rieu, 
1981). This is due to the fact that OH can be 
excited by infrared radiation from warm cir- 
cumstellar dust and that SiO is pumped by the 
stellar radiation. Strong SiS maser spikes also 
arise near the stellar surface. The pumping 
mechanism of H 2 0 masers is not well 
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Figure 3-3. The OH spectrum ( 1612-MHz satellite line) (top) and the corresponding 
maps at different radial velocities (with respect to the stellar velocity taken equal to 
zero) of the OH/IR source , OH26.5 + 0.6 , observed with the VLA by Baud (1981). 
The dot in each frame corresponds to the stellar position. Radial velocities in km 
s' 1 are noted on top (left) of each frame , and the beam shape is represented in the 
frame at bottom right. In a spherical shell, the maps corresponding to the extreme 
radial velocities , 15.1 and -16.7 km s' 1 (red-shifted and blue-shifted peaks , respec- 
tively), are associated with the part of the envelope confined within a small area in 
the double cone (see also Figure 3-2), where the maser amplification is maximum. 
The maps taken at velocities closer to the stellar velocity and corresponding to more 
extended areas outside the cone show a more complicated shape , in particular a ring 
structure (e.g., at velocities 3.8 and -5.3 km s' 1 ). 
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understood because of the complexity of the 
energy level structure. 

Figure 3-4 displays the SiS maser spectrum 
in the ground state (v = 0, J = 1-0 at 1.65-cm 
wavelength of the carbon star, IRC + 10216, 
observed at minimum light with the 140- ft radio 
telescope at Green Bank together with the 
calculated profile (dots) derived from an excita- 
tion model of the kind described in the section 
Multilevel Radiative Transfer (Nguyen-Q-Rieu 
et al., 1984). The narrow maser spike at -40 km 
s' 1 is likely to be produced by the amplifica- 
tion of the emission of the infrared core by the 
gas in the immediate vicinity of the core. The 
observations by Henkel et al. (1983) made at 
maximum infrared light indicate that the maser 
intensity is increased by - 40 percent as com- 
pared with the intensity at minimum. 

The increase of the infrared flux and of the 
mass loss during maximum light could be 
responsible for the enhancement of the SiS 
maser emission. It is noteworthy that weak CO 
maser emission can exist in the circumstellar 
envelope (Morris, 1980). 


Strong circumstellar masers are probably 
saturated. They are usually not subject to er- 
ratic time variation, but instead vary smoothly 
with the stellar phase. 

OH emission is known to be polarized, up 
to -60 percent (Robinson et al., 1970). 
Circular polarization is usually observed, but 
in some cases, polarization can change from cir- 
cular to linear. Significant linear polarization 
has been detected in SiO maser spectra (Troland 
et al., 1979). The position angles of the po- 
larized features vary with time. In the case of 
the Mira variable, R Cas, sudden changes in the 
position angles have been detected near the 
maximum of the optical light curve (Clark et 
al., 1982). Shock waves which can affect the 
gain paths and the orientation of the maser cells 
are believed to be responsible for this sudden 
variation. 

THERMAL EMISSION 

Thermal CO emission in the ground-state (u 
= 0, J = 1-0) at 2.6 mm has been discovered 



Figure 3-4. The 18155-MHz SiS (J = 1-0) spectrum of the carbon star IRC + 10216. Dots corre- 
spond to the theoretical spectrum derived from an excitation model. The blue-shifted feature at -40 
km s~ l is the result of an amplification of the background stellar core by the SiS circumstellar maser 
cloud (from Nguyen-Q-Rieu et al., 1984). 
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in the envelope of the carbon star, IRC + 
10216 (Solomon et al., 1971). Thermal emission 
from SiO, CO, CN, and CS (ground-state) and 
their isotopes was later detected in a number 
of oxygen-rich and carbon-rich stars, as well as 
in some S-type stars (Wilson et al., 1971; 
Zuckerman et al., 1977, 1978; Lo and Bechis, 
1977; Lambert and Vanden Bout, 1978; Wan- 
nier et al., 1979; Olofsson et al., 1982a). 
Molecular species like HNC, HCN, HC 3 N, 
and other linear carbon chain molecules also 
exist in the shell of carbon-rich stars (Jewell and 
Snyder, 1984; Olofsson et al., 1982b). 

Figure 3-5 shows the spectra of HNC (7 = 
1-0), SiS (7 = 5-4), and HC 3 N (7 = 10-9) 
observed simultaneously in the direction of IRC 
+ 10216 using the millimeter Onsala 

radiotelescope equipped with a broadband 
receiver (Olofsson et al., 1982a). 

Contrary to SiO masers which arise in the 
vicinity of the central star, thermal molecular 
line emission extends over the entire envelope. 
In this respect, thermal spectra give informa- 
tion on the physical parameters of the entire 
shell (namely, the mass loss, gas, and dust den- 
sities). For example, spatial CO mapping of the 
envelope of IRC + 10216 suggests that the 
present mass loss rate may have been higher in 



the past (Wannier et al., 1979). The deter- 
mination of the abundance of molecules and 
their isotopic substitutes is useful in the un- 
derstanding of the processing of matter involv- 
ing the mixing mechanisms (Iben, 1981) and the 
photodestruction process caused by the in- 
terstellar ultraviolet radiation field (Huggins 
and Glassgold, 1982). Neutral/neutral reactions 
seem to be important in circumstellar shells, 
while interstellar chemistry is dominated by 
ion/molecule reactions. This is due to the fact 
that high temperatures ( >100 K) allow endo- 
thermic reactions, as well as others with activa- 
tion energy barriers, to proceed in the circum- 
stellar envelope (Scalo and Slavsky, 1980; La- 
font et al., 1982). 

Molecular Excitation and 
Line Profiles 

The excitation of circumstellar molecules is 
governed by Equation (3-8) in which both col- 
lisional and radiative processes are involved. 
Collisions, essentially with H 2 , tend to ther- 
malize the lines. This is usually the case with 
CO lines in dense envelopes. Infrared radiation 
from the central star and dust embedded in the 
envelope can produce population inversion as 
discussed in the section Characteristics of Cir- 
cumstellar Masers . 

The main-beam brightness temperature, 
T mb , at a radial velocity, K , is obtained by 
convolving the brightness temperature distribu- 
tion, T b (see Equation (3-15)) with the an- 
tenna main beam, which is assumed to be gaus- 
sian (Olofsson et al., 1982a): 


T M„ (*o = 


81n2 


MB v > g2 


cm ax 

L 


pmin ( 3 - 16 ) 


T R (p, K)e' 41n2 (P/B) pdp , 


Figure 3-5 . Simultaneous observations of HNC, 
SiS, and HC yV lines around 91000 MHz from 
IRC 4- 10216, using a broadband receiver in- 
stalled at the focus of the 20-m Onsala radio 
telescope (from Olofsson et al., 1982a). 


where B is the full width at half power (FWHP) 
of the antenna beam and p = r j 1 - 
(V/V ex p) 2 \ 1/2 (see Figure 3-1). In Equation 
(3-16), we have assumed a uniformly excited 
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spherical shell (T is constant throughout the 
shell), expanding at a constant velocity. If the 
mass-loss rate is constant, the total particle den- 
sity is N(r) oc f 2 . The optical depth can be ex- 
pressed as (Olofsson et al., 1982a): 

t(P,V) = r(B,0) X B 

{p 

The diagram (Figure 3-6) which represents the 
line profiles T MB (V r ) as a function of the op- 
tical depth, t (B,0), and the shell size, R e , ex- 
pressed in units of beamwidth, B, is calculated 
using Equation (3-16). 

For an unresolved optically thin shell (posi- 
tion 1 in the diagram), the profile is rectangular. 
The line shape is parabolic if the shell is optical- 
ly thick (position 2) and becomes flatter when 
the shell size is extended with respect to the 



antenna beam (positions 3 and 5). The profile 
can exhibit a double-peaked structure for an 
optically thin and resolved shell (position 4). 
These effects can be explained by the fact that, 
when the source is resolved by the telescope 
beam, the part of the shell which is not inside 
the beam corresponds to low radial velocity line 
center. As a result, the contribution of this 
outer shell to the main-beam brightness temper- 
ature is not important and leads to a more or 
less pronounced depression at the line center. 
Therefore, a quick inspection of the line pro- 
files can indicate the nature of the circumstellar 
shell. 

Line Asymmetry and Time Variation 

Figure 3-7 shows a variety of line profiles 
observed in the direction of IRC + 10216 



R e /B 


Figure 3-6. A diagram representing various line shapes as a function of optical depth 
(ordinate) and source size R e in unit of beamwidth B (abscissa). Dots correspond to 
the points in the diagram for which the line profiles are calculated. 


I 
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Figure 3-7. A sample of spectra (CO at 115271 MHz , HNC at 90663 MHz , and HC 3 N at 90979 MHz) 
from IRC + 10216 obtained with a high resolution spectrometer v/tsv ~ 4 x 10 5 ) (from Olofsson 
et al. f 1982a). V rCD is the radial velocity relative to the local standard of rest. 

LJA 


(Olofsson et al., 1982a). The spectra correspond 
to different line opacities and source sizes, ac- 
cording to the diagram shown in Figure 3-6. 
The CO line emanates from an optically thick 
extended shell. The HNC and HC 3 N transi- 
tions correspond to optically thin regions. 
Although the HNC shell is resolved by the 
telescope beam, the HC 3 N shell appears to be 
a point source. Theoretically, an estimate of the 
shell parameters can be derived by fitting the 
main-beam brightness temperature profiles, 
calculated by using Equation (3-16), to the 
observed spectra. However, the observed pro- 
files (Figure 3-7) present some asymmetry in the 
sense that the blue-shifted wing is systematically 
weaker than the red-shifted one. This effect is 
most prominent in the case of 13 CO (Figure 
3-8). It should be mentioned that in the SiS J 


= 1-0 profile (Figure 3-4), the depression of 
the blue-shifted wing is disguised by the maser 
action which amplifies the stellar radiation and 
thereby gives rise to a narrow spike at the blue 
wing. 

In fact, the hypothesis of a constant excita- 
tion temperature throughout the circumstellar 
cloud cannot account for the line asymmetry. 
The inner shell is more excited and hotter than 
the outer shell. As a result, along the line of 
sight, outer colder material is absorbing emis- 
sion from the hotter inner layers in the front 
hemisphere, but such absorption does not af- 
fect the back hemisphere. Similar asymmetry 
also seems to be observed in a Ori (Knapp et 
al., 1980). 

In order to account for this effect, instead 
of using Equation (3-15), we shall calculate the 
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Figure 3-8. The 13 CO spectrum (110201 MHz) 
from IRC + 10216 showing the line asymmetry 
(blue-shifted wing weaker than red-shifted 
wing). Dots represent the spectrum derived 
from a theoretical model (Nguyen-Q-Rieu et 
al. y 1984) which explains the asymmetry in 
terms of an absorption of the emission from 
the hot inner part of the envelope by the colder 
outer layers. 


brightness temperature separately for the front 
(T b ) and back (T B+ ) hemispheres by in- 
tegrating along the line of sight (see Figure 3-1): 


k(z # ) dz # , 

x(0 is the line absorption coefficient. 

Figure 3-8 displays the observed 13 CO spec- 
trum of IRC + 10216, together with the spec- 
trum computed (dots) by Equations (3-17), 
(3-18), and (3-16) (Nguyen-Q-Rieu et al., 1984; 
see also Morris et al., 1985). 

The molecular line shape is sensitive to the 
variation of infrared light. At the maximum of 
the infrared light curve, infrared photons ( - 10 
/im) escaping from the central heating sources 
(star or infrared core) can affect the outer shell 
and enhance the line wings (Sahai et al., 1984). 
The contrast between the line wings and the line 
center can be higher at maximum light than at 
minimum. In the case of IRC + 10216, it has 
been shown that, at maximum, the line center 
and the line wings of the SiS J = 6-5 transi- 
tion can be enhanced by - 20 and - 40 percent, 
respectively (Nguyen-Q-Rieu et al., 1984). 

CONCLUSION AND PROSPECTS 



T 

1 B- 


(3-17) 


+S — oo 



Radio molecular lines appear to be useful 
probes into the stellar environment. SiO masers 
which are excited by stellar photons to a vibra- 
tional level as high as v = 3 (Scalise and Lepine, 
1978), which corresponds to an energy E/k - 
5300 K, provide information on the physical 
conditions in the immediate vicinity of the 
stellar photosphere. Lower vibrationally excited 
transitions of SiO occur in u = 0 and 1 (E/k 
< 1 800 K) sample regions in which dust grains 
begin to form. The proximity of SiO masers to 
the pulsating stellar atmosphere implies that 
SiO observations can shed light on the kine- 
matics and mass-loss process in the innermost 
layers. As explained in the sections Multilevel 
Radiative Transfer and Characteristics of Cir- 
cumstellar Masers, the maser amplification and 
hence the line shape depend on the velocity gra- 
dient. SiO maser emission does not usually ex- 
hibit a double-peaked profile typical of OH 
masers, but a single feature centered at the 
stellar velocity. This observational difference 
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suggests that the velocity gradient in the inner 
region in which SiO occurs is higher (see Equa- 
tion (3-13), where y = d\og V/dlogr is now 
>> 1) than in the external envelope in which 
OH masers are located (7 « 1). OH masers 
and thermal molecular emission of CO and 
other molecules, including cyanopolyynes, 
which can be excited by collisions with H 2 and 
by infrared radiation from circumstellar dust, 
can be used as probes into the physical condi- 
tions (on a larger scale) throughout the 
envelope. Valuable information on the physics 
operating in the envelope of IRC + 10216 has 
been recently obtained by high sensitivity obser- 
vations and detailed theoretical analyses (Mor- 
ris, 1975; Wannier et al., 1979; Kwan and 
Linke, 1982; Olofsson, et al., 1982a; Nguyen- 
Q-Rieu et al., 1984). However, accurate cir- 
cumstellar chemistry still suffers from the lack 
of data on chemical processes. Chemistry is 
usually assumed to be frozen at its equilibrium 
state in the inner region and unchanging dur- 
ing the expansion of the envelope (McCabe et 
al., 1979). This is probably true for species 
which are not highly reactive, but not for 
radicals. Nonequilibrium chemistry may take 
place, for instance, in the regions in which 
molecular abundances are governed by photo- 
dissociation, which tends to produce radicals 
(Goldreich and Scoville, 1976; Scalo and Slav- 
sky, 1980; Lafont et al., 1982). In this respect, 
mapping of the envelope of IRC + 10216 in 
the molecular lines emitted by some radicals 
and parent molecules is highly desirable. Pre- 
sent molecular observations of the envelopes of 
stars of various types (Olofsson et al., 1982b) 
and those with the millimeter radio telescopes 
of Institut de Radioastronomie Millimetrique 
(IRAM) may also cast some light on the prob- 
lem of chemistry (see Glassgold and Huggins, 
this volume). 

Infrared speckle interferometry in the 
molecular lines (Dyck et al., 1983) and in the 
continuum (Foy et al., 1979; Dyck et al., 1984) 
is helpful in the investigation of the inner region 
of the envelope. 

The distribution of the circumstellar dust, 
whose emission contributes to the molecular 


excitation, can be determined by photometric 
measurements with a range of 1 to - 100 /un. 
In this respect, the preliminary results obtained 
with the instruments on board the Infrared 
Astronomical Satellite (IRAS) seem to indicate 
that many infrared objects associated with OH 
masers are extremely red and may experience 
mass-loss rates higher than 10“ 4 M© yr' 1 
(Olnon et al., 1984). High resolution spectro- 
scopy using the Infrared Space Observatory 
(ISO) will also be of great importance in the 
detection of molecular transitions in a wide 
range of infrared wavelengths. 
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